Autoantigen profiling reveals a shared post-COVID signature in fully recovered and Long COVID patients (preprint)

Some individuals do not return to baseline health following SARS-CoV-2 infection, leading to a condition known as Long COVID. The underlying pathophysiology of Long COVID remains unknown. Given that autoantibodies have been found to play a role in severity of COVID infection and certain other post-COVID sequelae, their potential role in Long COVID is important to investigate. Here we apply a well-established, unbiased, proteome-wide autoantibody detection technology (PhIP-Seq) to a robustly phenotyped cohort of 121 individuals with Long COVID, 64 individuals with prior COVID-19 who reported full recovery, and 57 pre-COVID controls. While a distinct autoreactive signature was detected which separates individuals with prior COVID infection from those never exposed to COVID, we did not detect patterns of autoreactivity that separate individuals with Long COVID relative to individuals fully recovered from SARS-CoV-2 infection. These data suggest that there are robust alterations in autoreactive antibody profiles due to infection; however, no association of autoreactive antibodies and Long COVID was apparent by this assay.

Systematic functional interrogation of SARS-CoV-2 host factors using Perturb-seq (preprint)

Numerous host factors of SARS-CoV-2 have been identified by screening approaches, but delineating their molecular roles during infection and whether they can be targeted for antiviral intervention remains a challenge. Here we use Perturb-seq, a single-cell CRISPR screening approach, to investigate how CRISPR interference of host factors changes the course of SARS-CoV-2 infection and the host response in human lung epithelial cells. Our data reveal two classes of host factors with pronounced phenotypes: factors required for the response to interferon and factors required for entry or early infection. Among the latter, we have characterized the NF-{kappa}B inhibitor I{kappa}B (NFKBIA), as well as the translation factors EIF4E2 and EIF4H as strong host dependency factors acting early in infection. Overall, our study provides high-throughput functional validation of host factors of SARS-CoV-2 and describes their roles during viral infection in both infected and bystander cells.

Autoantibody discovery across monogenic, acquired, and COVID19-associated autoimmunity with scalable PhIP-Seq (preprint)

Phage Immunoprecipitation-Sequencing (PhIP-Seq) allows for unbiased, proteome-wide autoantibody discovery across a variety of disease settings, with identification of disease-specific autoantigens providing new insight into previously poorly understood forms of immune dysregulation. Despite several successful implementations of PhIP-Seq for autoantigen discovery, including our previous work (Vazquez et al. 2020), current protocols are inherently difficult to scale to accommodate large cohorts of cases and importantly, healthy controls. Here, we develop and validate a high throughput extension of PhIP-seq in various etiologies of autoimmune and inflammatory diseases, including APS1, IPEX, RAG1/2 deficiency, Kawasaki Disease (KD), Multisystem Inflammatory Syndrome in Children (MIS-C), and finally, mild and severe forms of COVID19. We demonstrate that these scaled datasets enable machine-learning approaches that result in robust prediction of disease status, as well as the ability to detect both known and novel autoantigens, such as PDYN in APS1 patients, and intestinally expressed proteins BEST4 and BTNL8 in IPEX patients. Remarkably, BEST4 antibodies were also found in 2 patients with RAG1/2 deficiency, one of whom had very early onset IBD. Scaled PhIP-Seq examination of both MIS-C and KD demonstrated rare, overlapping antigens, including CGNL1, as well as several strongly enriched putative pneumonia-associated antigens in severe COVID19, including the endosomal protein EEA1. Together, scaled PhIP-Seq provides a valuable tool for broadly assessing both rare and common autoantigen overlap between autoimmune diseases of varying origins and etiologies.

Full genome Nobecovirus sequences from Malagasy fruit bats define a unique evolutionary history for this coronavirus clade (preprint)

Bats are natural reservoirs for both Alpha- and Betacoronaviruses and the hypothesized original hosts of five of seven known zoonotic coronaviruses. To date, the vast majority of bat coronavirus research has been concentrated in Asia, though coronaviruses are globally distributed; indeed, SARS-CoV and SARS-CoV-2-related Betacoronaviruses in the subgenus Sarbecovirus have been identified circulating in Rhinolophid bats in both Africa and Europe, despite the relative dearth of surveillance in these regions. As part of a long-term study examining the dynamics of potentially zoonotic viruses in three species of endemic Madagascar fruit bat (Pteropus rufus, Eidolon dupreanum, Rousettus madagascariensis), we carried out metagenomic Next Generation Sequencing (mNGS) on urine, throat, and fecal samples obtained from wild-caught individuals. We report detection of RNA derived from Betacoronavirus subgenus Nobecovirus in fecal samples from all three species and describe full genome sequences of novel Nobecoviruses in P. rufus and R. madagascariensis. Phylogenetic analysis indicates the existence of five distinct Nobecovirus clades, one of which is defined by the highly divergent sequence reported here from P. rufus bats. Madagascar Nobecoviruses derived from P. rufus and R. madagascariensis demonstrate, respectively, Asian and African phylogeographic origins, mirroring those of their fruit bat hosts. Bootscan recombination analysis indicates significant selection has taken place in the spike, nucleocapsid, and NS7 accessory protein regions of the genome for viruses derived from both bat hosts. Madagascar offers a unique phylogeographic nexus of bats and viruses with both Asian and African phylogeographic origins, providing opportunities for unprecedented mixing of viral groups and, potentially, recombination. As fruit bats are handled and consumed widely across Madagascar for subsistence, understanding the landscape of potentially zoonotic coronavirus circulation is essential for mitigation of future zoonotic threats.

Using genomic epidemiology of SARS-CoV-2 to support contact tracing and public health surveillance in rural Humboldt County, California. (preprint)

During the COVID-19 pandemic within the United States, much of the responsibility for diagnostic testing and epidemiologic response has relied on the action of county-level departments of public health. Here we describe the integration of genomic surveillance into epidemiologic response within Humboldt County, a rural county in northwest California. Through a collaborative effort, 853 whole SARS-CoV-2 genomes were generated, representing [~]58% of the 1,449 SARS-CoV-2-positive cases detected in Humboldt County as of mid-March 2021. Phylogenetic analysis of these data was used to develop a comprehensive understanding of SARS-CoV-2 introductions to the county and to support contact tracing and epidemiologic investigations of all large outbreaks in the county. In the case of an outbreak on a commercial farm, viral genomic data were used to validate reported epidemiologic links and link additional cases within the community who did not report a farm exposure to the outbreak. During a separate outbreak within a skilled nursing facility, genomic surveillance data were used to rule out the putative index case, detect the emergence of an independent Spike:N501Y substitution, and verify that the outbreak had been brought under control. These use cases demonstrate how developing genomic surveillance capacity within local public health departments can support timely and responsive deployment of genomic epidemiology for surveillance and outbreak response based on local needs and priorities.

Upper airway gene expression reveals a more robust innate and adaptive immune response to SARS-CoV-2 in children compared with older adults (preprint)

Unlike other respiratory viruses, SARS-CoV-2 disproportionately causes severe disease in older adults and only rarely in children. To investigate whether differences in the upper airway immune response could contribute to this disparity, we compared nasopharyngeal gene expression in 83 children (<19-years-old; 38 with SARS-CoV-2, 11 with other respiratory viruses, 34 with no virus) and 154 adults (>40-years-old; 45 with SARS-CoV-2, 28 with other respiratory viruses, 81 with no virus). Expression of interferon-stimulated genes (ISGs) was robustly activated in both children and adults with SARS-CoV-2 compared to the respective non-viral groups, with only relatively subtle distinctions. Children, however, demonstrated markedly greater upregulation of pathways related to B cell and T cell activation and proinflammatory cytokine signaling, including TNF, IFN{gamma}, IL-2 and IL-4 production. Cell type deconvolution confirmed greater recruitment of B cells, and to a lesser degree macrophages, to the upper airway of children. Only children exhibited a decrease in proportions of ciliated cells, the primary target for SARS-CoV-2, upon infection with the virus. These findings demonstrate that children elicit a more robust innate and adaptive immune response to SARS-CoV-2 infection in the upper airway that likely contributes to their protection from severe disease in the lower airway.

Impaired antibacterial immune signaling and changes in the lung microbiome precede secondary bacterial pneumonia in COVID-19 (preprint)

Secondary bacterial infections, including ventilator associated pneumonia (VAP), lead to worse clinical outcomes and increased mortality following viral respiratory infections. Critically ill patients with coronavirus disease 2019 (COVID-19) face an elevated risk of VAP, although susceptibility varies widely. Because mechanisms underlying VAP predisposition remained unknown, we assessed lower respiratory tract host immune responses and microbiome dynamics in 36 patients, including 28 COVID-19 patients, 15 of whom developed VAP, and eight critically ill controls. We employed a combination of tracheal aspirate bulk and single cell RNA sequencing (scRNA-seq). Two days before VAP onset, a lower respiratory transcriptional signature of bacterial infection was observed, characterized by increased expression of neutrophil degranulation, toll-like receptor and cytokine signaling pathways. When assessed at an earlier time point following endotracheal intubation, more than two weeks prior to VAP onset, we observed a striking early impairment in antibacterial innate and adaptive immune signaling that markedly differed from COVID-19 patients who did not develop VAP. scRNA-seq further demonstrated suppressed immune signaling across monocytes/macrophages, neutrophils and T cells. While viral load did not differ at an early post-intubation timepoint, impaired SARS-CoV-2 clearance and persistent interferon signaling characterized the patients who later developed VAP. Longitudinal metatranscriptomic analysis revealed disruption of lung microbiome community composition in patients who developed VAP, providing a connection between dysregulated immune signaling and outgrowth of opportunistic pathogens. Together, these findings demonstrate that COVID-19 patients who develop VAP have impaired antibacterial immune defense weeks before secondary infection onset.

Longitudinal single-cell epitope and RNA-sequencing reveals the immunological impact of type 1 interferon autoantibodies in critical COVID-19 (preprint)

Type I interferon (IFN-I) neutralizing autoantibodies have been found in some critical COVID-19 patients; however, their prevalence and longitudinal dynamics across the disease severity scale, and functional effects on circulating leukocytes remain unknown. Here, in 284 COVID-19 patients, we found IFN-I autoantibodies in 19% of critical, 6% of severe and none of the moderate cases. Longitudinal profiling of over 600,000 peripheral blood mononuclear cells using multiplexed single-cell epitope and transcriptome sequencing from 54 COVID-19 patients, 15 non-COVID-19 patients and 11 non-hospitalized healthy controls, revealed a lack of IFN-I stimulated gene (ISG-I) response in myeloid cells from critical cases, including those producing anti-IFN-I autoantibodies. Moreover, surface protein analysis showed an inverse correlation of the inhibitory receptor LAIR-1 with ISG-I expression response early in the disease course. This aberrant ISG-I response in critical patients with and without IFN-I autoantibodies, supports a unifying model for disease pathogenesis involving ISG-I suppression via convergent mechanisms.

Estimation of secondary household attack rates for emergent SARS-CoV-2 variants detected by genomic surveillance at a community-based testing site in San Francisco (preprint)

BackgroundSequencing of the SARS-CoV-2 viral genome from patient samples is an important epidemiological tool for monitoring and responding to the pandemic, including the emergence of new mutations in specific communities. MethodsSARS-CoV-2 genomic sequences were generated from positive samples collected, along with epidemiological metadata, at a walk-up, rapid testing site in the Mission District of San Francisco, California during November 22-December 2, 2020 and January 10-29, 2021. Secondary household attack rates and mean sample viral load were estimated and compared across observed variants. ResultsA total of 12,124 tests were performed yielding 1,099 positives. From these, 811 high quality genomes were generated. Certain viral lineages bearing spike mutations, defined in part by L452R, S13I, and W152C, comprised 54.9% of the total sequences from January, compared to 15.7% in November. Household contacts exposed to "West Coast" variants were at higher risk of infection compared to household contacts exposed to lineages lacking these variants (0.357 vs 0.294, RR=1.29; 95% CI:1.01-1.64). The reproductive number was estimated to be modestly higher than other lineages spreading in California during the second half of 2020. Viral loads were similar among persons infected with West Coast versus non-West Coast strains, as was the proportion of individuals with symptoms (60.9% vs 64.1%). ConclusionsThe increase in prevalence, relative household attack rates, and reproductive number are consistent with a modest transmissibility increase of the West Coast variants; however, additional laboratory and epidemiological studies are required to better understand differences between these variants. SummaryWe observed a growing prevalence and elevated attack rate for "West Coast" SARS-CoV-2 variants in a community testing setting in San Francisco during January 2021, suggesting its modestly higher transmissibility.

COVID-19 ARDS is characterized by a dysregulated host response that differs from cytokine storm and may be modified by dexamethasone (preprint)

We performed comparative lower respiratory tract transcriptional profiling of 52 critically ill patients with ARDS from COVID-19 or other etiologies, or without ARDS. We found no evidence of cytokine storm but instead observed complex host response dysregulation driven by genes with non-canonical roles in inflammation and immunity that were predicted to be modulated by dexamethasone. Compared to other viral ARDS, COVID-19 was characterized by impaired interferon-stimulated gene expression.

Direct detection of SARS-CoV-2 using CRISPR-Cas13a and a mobile phone (preprint)

The December 2019 outbreak of a novel respiratory virus, SARS-CoV-2, has become an ongoing global pandemic due in part to the challenge of identifying symptomatic, asymptomatic and pre- symptomatic carriers of the virus. CRISPR-based diagnostics that utilize RNA and DNA-targeting enzymes can augment gold-standard PCR-based testing if they can be made rapid, portable and accurate. Here we report the development of an amplification-free CRISPR-Cas13a-based mobile phone assay for direct detection of SARS-CoV-2 from nasal swab RNA extracts. The assay achieved ~100 copies/L sensitivity in under 30 minutes and accurately detected a set of positive clinical samples in under 5 minutes. We combined crRNAs targeting SARS-CoV-2 RNA to improve sensitivity and specificity, and we directly quantified viral load using enzyme kinetics. Combined with mobile phone-based quantification, this assay can provide rapid, low-cost, point-of-care screening to aid in the control of SARS-CoV-2.

Nasopharyngeal SARS-CoV2 viral loads in young children do not differ significantly from those in older children and adults (preprint)

The role of children in the spread of the SARS-CoV2 coronavirus has become a matter of urgent debate as societies in the US and abroad consider how to safely reopen schools. Small studies have suggested higher viral loads in young children. Here we present a multicenter investigation on over five thousand SARS-CoV-2 cases confirmed by real-time reverse transcription (RT) PCR assay. Notably, we found no discernable difference in amount of viral nucleic acid among young children and adults.

Immunologically distinct responses occur in the CNS of COVID-19 patients (preprint)

One third of COVID-19 patients develop significant neurological symptoms, yet SARS-CoV-2 is rarely detected in central nervous system (CNS) tissue, suggesting a potential role for parainfectious processes, including neuroimmune responses. We therefore examined immune parameters in cerebrospinal fluid (CSF) and blood samples from a cohort of patients with COVID-19 and significant neurological complications. We found divergent immunological responses in the CNS compartment, including increased levels of IL-12 and IL-12-associated innate and adaptive immune cell activation. Moreover, we found increased proportions of B cells in the CSF relative to the periphery and evidence of clonal expansion of CSF B cells, suggesting a divergent intrathecal humoral response to SARS-CoV-2. Indeed, all COVID-19 cases examined had anti-SARS-CoV-2 IgG antibodies in the CSF whose target epitopes diverged from serum antibodies. We directly examined whether CSF resident antibodies target self-antigens and found a significant burden of CNS autoimmunity, with the CSF from most patients recognizing neural self-antigens. Finally, we produced a panel of monoclonal antibodies from patients CSF and show that these target both anti-viral and anti-neural antigens--including one mAb specific for the spike protein that also recognizes neural tissue. This exploratory immune survey reveals evidence of a compartmentalized and self-reactive immune response in the CNS meriting a more systematic evaluation of neurologically impaired COVID-19 patients. One Sentence SummaryA subset of COVID-19 patients with neurologic impairment show cerebrospinal fluid-specific immune alterations that point to both neuroinvasion and anti-neural autoimmunity as potential causes of impairment.

Identification of a polymorphism in the N gene of SARS-CoV-2 that adversely impacts detection by a widely-used RT-PCR assay (preprint)

We identify a mutation in the N gene of SARS-CoV-2 that adversely affects annealing of a commonly used RT-PCR primer; epidemiologic evidence suggests the virus retains pathogenicity and competence for spread. This reinforces the importance of using multiple targets, preferably in at least 2 genes, for robust SARS-CoV-2 detection. Article Summary LineA SARS-CoV-2 variant that occurs worldwide and has spread in California significantly affects diagnostic sensitivity of an N gene assay, highlighting the need to employ multiple viral targets for detection.

The SARS-CoV-2 Spike mutation D614G increases entry fitness across a range of ACE2 levels, directly outcompetes the wild type, and is preferentially incorporated into trimers (preprint)

Early in the current pandemic, the D614G mutation arose in the Spike protein of SARS-CoV-2 and quickly became the dominant variant globally. Mounting evidence suggests D614G enhances viral entry. Here we use a direct competition assay with single-cycle viruses to show that D614G outcompetes the wildtype. We developed a cell line with inducible ACE2 expression to confirm that D614G more efficiently enters cells with ACE2 levels spanning the different primary cells targeted by SARS-CoV-2. Using a new assay for crosslinking and directly extracting Spike trimers from the pseudovirus surface, we found an increase in trimerization efficiency and viral incorporation of D614G protomers. Our findings suggest that D614G increases infection of cells expressing a wide range of ACE2, and informs the mechanism underlying enhanced entry. The tools developed here can be broadly applied to study other Spike variants and SARS-CoV-2 entry, to inform functional studies of viral evolution and vaccine development.

Furin Cleavage Site Is Key to SARS-CoV-2 Pathogenesis (preprint)

SARS-CoV-2 has resulted in a global pandemic and shutdown economies around the world. Sequence analysis indicates that the novel coronavirus (CoV) has an insertion of a furin cleavage site (PRRAR) in its spike protein. Absent in other group 2B CoVs, the insertion may be a key factor in the replication and virulence of SARS-CoV-2. To explore this question, we generated a SARS-CoV-2 mutant lacking the furin cleavage site ({Delta}PRRA) in the spike protein. This mutant virus replicated with faster kinetics and improved fitness in Vero E6 cells. The mutant virus also had reduced spike protein processing as compared to wild-type SARS-CoV-2. In contrast, the {Delta}PRRA had reduced replication in Calu3 cells, a human respiratory cell line, and had attenuated disease in a hamster pathogenesis model. Despite the reduced disease, the {Delta}PRRA mutant offered robust protection from SARS-CoV-2 rechallenge. Importantly, plaque reduction neutralization tests (PRNT50) with COVID-19 patient sera and monoclonal antibodies against the receptor-binding domain found a shift, with the mutant virus resulting in consistently reduced PRNT50 titers. Together, these results demonstrate a critical role for the furin cleavage site insertion in SARS-CoV-2 replication and pathogenesis. In addition, these findings illustrate the importance of this insertion in evaluating neutralization and other downstream SARS-CoV-2 assays. ImportanceAs COVID-19 has impacted the world, understanding how SARS-CoV-2 replicates and causes virulence offers potential pathways to disrupt its disease. By removing the furin cleavage site, we demonstrate the importance of this insertion to SARS-CoV-2 replication and pathogenesis. In addition, the findings with Vero cells indicate the likelihood of cell culture adaptations in virus stocks that can influence reagent generation and interpretation of a wide range of data including neutralization and drug efficacy. Overall, our work highlights the importance of this key motif in SARS-CoV-2 infection and pathogenesis. Article SummaryA deletion of the furin cleavage site in SARS-CoV-2 amplifies replication in Vero cells, but attenuates replication in respiratory cells and pathogenesis in vivo. Loss of the furin site also reduces susceptibility to neutralization in vitro.

Epitopes targeted by T cells in convalescent COVID-19 patients (preprint)

There is growing evidence pointing to the protective role of T cells against COVID-19. Vaccines eliciting targeted T cell responses have the potential to provide robust, long-lasting immunity. However, their design requires knowledge of the SARS-CoV-2-specific epitopes that can elicit a T cell response and confer protection across a wide population. Here, we provide a unified description of emerging data of SARS-CoV-2 T cell epitopes compiled from results of 8 independent studies of convalescent COVID-19 patients. We describe features of these epitopes relevant for vaccine design, while indicating knowledge gaps that can, in part, be augmented using prior immunological data from SARS-CoV. The landscape of SARS-CoV-2 T cell epitopes that we describe can help guide SARS-CoV-2 vaccine development as well as future immunological studies. A web-based platform has also been developed to complement these efforts.

AI aided design of epitope-based vaccine for the induction of cellular immune responses against SARS-CoV-2 (preprint)

The heavy burden imposed by the COVID-19 pandemic on our society triggered the race towards the development of therapies or preventive strategies. Among these, antibodies and vaccines are particularly attractive because of their high specificity, low probability of drug-drug interaction, and potentially long-standing protective effects. While the threat at hand justifies the pace of research, the implementation of therapeutic strategies cannot be exempted from safety considerations. There are several potential adverse events reported after the vaccination or antibody therapy, but two are of utmost importance: antibody-dependent enhancement (ADE) and cytokine storm syndrome (CSS). On the other hand, the depletion or exhaustion of T-cells has been reported to be associated with worse prognosis in COVID-19 patients. This observation suggests a potential role of vaccines eliciting cellular immunity, which might simultaneously limit the risk of ADE and CSS. Such risk was proposed to be associated with FcR-induced activation of proinflammatory macrophages (M1) by Fu et al. 2020 and Iwasaki et al. 2020. All aspects of the newly developed vaccine (including the route of administration, delivery system, and adjuvant selection) may affect its effectiveness and safety. In this work we use a novel in silico approach (based on AI and bioinformatics methods) developed to support the design of epitope-based vaccines. We evaluated the capabilities of our method for predicting the immunogenicity of epitopes. Next, the results of our approach were compared with other vaccine-design strategies reported in the literature. The risk of immuno-toxicity was also assessed. The analysis of epitope conservation among other Coronaviridae was carried out in order to facilitate the selection of peptides shared across different SARS-CoV-2 strains and which might be conserved in emerging zootic coronavirus strains. Finally, the potential applicability of the selected epitopes for the development of a vaccine eliciting cellular immunity for COVID-19 was discussed, highlighting the benefits and challenges of such an approach.

Upper airway gene expression differentiates COVID-19 from other acute respiratory illnesses and reveals suppression of innate immune responses by SARS-CoV-2 (preprint)

We studied the host transcriptional response to SARS-CoV-2 by performing metagenomic sequencing of upper airway samples in 238 patients with COVID-19, other viral or non-viral acute respiratory illnesses (ARIs). Compared to other viral ARIs, COVID-19 was characterized by a diminished innate immune response, with reduced expression of genes involved in toll-like receptor and interleukin signaling, chemokine binding, neutrophil degranulation and interactions with lymphoid cells. Patients with COVID-19 also exhibited significantly reduced proportions of neutrophils and macrophages, and increased proportions of goblet, dendritic and B-cells, compared to other viral ARIs. Using machine learning, we built 26-, 10- and 3-gene classifiers that differentiated COVID-19 from other acute respiratory illnesses with AUCs of 0.980, 0.950 and 0.871, respectively. Classifier performance was stable at low viral loads, suggesting utility in settings where direct detection of viral nucleic acid may be unsuccessful. Taken together, our results illuminate unique aspects of the host transcriptional response to SARS-CoV-2 in comparison to other respiratory viruses and demonstrate the feasibility of COVID-19 diagnostics based on patient gene expression.

ReScan, a Multiplex Diagnostic Pipeline, Pans Human Sera for SARS-CoV-2 Antigens (preprint)

Serologic assays are needed to determine SARS-CoV-2 seroprevalence, but poor specificity can overestimate exposures. Here, we built a pan-human coronavirus proteome-wide programmable phage display assay (VirScan) to profile coronavirus antigens specifically enriched by 20 COVID-19 patient serum IgG. With ReScan, a new diagnostic development workflow which combines the isolation of phage expressing the most immunogenic peptides with paper-based microarrays manufactured via acoustic liquid handling, we identified 9 candidate antigens from a library of 534 SARS-CoV-2 peptides. These arrays could form the basis of a multiplexed COVID-19 serologic assay with enhanced specificity. ReScan has broad applicability for serologic assay development.